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Abstract 

We analyse in detail the possibility of Bose-Einstein condensation of pions produced in heavy-ion 
collisions at the beam energy ^snn = 2.76 TeV. Our approach is based on the chemical non-equilibrium 
thermal model of hadron production which has been generalised to include separately the contribution 
from the local zero-momentum state. In order to study both the hadronic multiplicities and the 
transverse-momentum spectra, we use the Cracow freeze-out model which parameterises the flow and 
space-time geometry of the system at freeze-out in a very economic way. Our analysis indicates that 
about 5% of all pions may form the Bose-Einstein condensate. 
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I. INTRODUCTION 


Thermal models of hadron production serve us as the basic tool to obtain general information 
about the properties of matter produced in heavy-ion collisions [IH25]. The successes of such 
models in description of different collision processes indicate typically that the produced matter 
is well thermalised and exhibits collective behavior. This, in turn, suggests the formation of 
the equilibrated system of quarks and gluons at earlier stages of the collisions. 

The newest LHC data on hadron production indicates that the most common thermal ap¬ 
proach, based on the chemical and thermal equilibration assumption, fails to describe the ratio 
of protons to pions [26] . The mean multiplicities of protons and anti-protons also deviate from 
the £t for about three standard deviations [27] . In addition, most of the hydrodynamic calcula¬ 
tions, that have been successfully used to reproduce the harmonic flow coefficients, have prob¬ 
lems with the correct predictions for very low transverse-momentum spectra of pions [221 12B1 - I5U] . 

In our previous works [SUES] we have demonstrated that these two problems may be solved 
by the assumption that matter produced in heavy-ion collisions at the LHC energies is formed 
out of chemical equilibrium min]. One of the predictions of the chemical non-equilibrium 
model is that the pion abundances are characterized by the non-zero value of the chemical 
potentiah which is very close to the critical value for the Bose-Einstein condensation (BEC). 
This suggests a possible onset of BEC at the LHC energies. 

The recent analysis of two- and three-pion correlations done by ALICE collaboration [M] 
shows an intriguing result that a coherent fraction in charged pion emission may reach 23%. 
The main aim of our present work is to determine quantitatively the amount of particles in 
condensate at different centralities using mean multiplicities and spectra of the particles. 

The fact that the pion chemical potential is very close to its critical value requires that 
thermal analyses in this case should be performed in a more careful way with the explicit 
treatment of the local zero-momentum state [35]. As the inclusion of the pion ground state 
effects is relatively easy in the analyses of the ratios of hadronic abundances, a fully covariant 
framework describing local formation of the condensate in an expanding medium is not at hand. 
The formulation of such a framework requires that the Cooper-Frye formula is generalized in 
a very special way. In this work we demonstrate how this can be achieved for boost-invariant 
and cylindrically symmetric systems. We make the corresponding changes in the previously 

^ The equilibrium model [53] also shows the hints of pion chemical potential at the LHC. 
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used by us SHARE [SB] and THERMINATOR [371 |3B] codes. 

The hadronic data are analyzed within three possible scenarios: using the full equilibrium 
model (EQ), the chemical non-equilibrium model (NEQ) that have been used in [311132], and 
the newly developed chemical non-equilibrium model with Bose-Einstein condensate (BEC) 
We include all available LHC data on Pb-|-Pb collisions at ^snn = 2.76 TeV [2Sl 
in our study using the modihed SHARE and use the modihed THERMINATOR in order to 
check the effect of BEC on the pion spectra at low transverse-momenta [261 112]. We also 
discuss in greater detail the effects which the data on neutral pion production may have on 
our results. Our comparisons between different thermal frameworks indicate the plausibility of 
pion condensation in heavy-ion collisions at the LHC. 

The paper is organized as follows: in Section [^ we describe the role of the ground state and 


the difference between EQ, NEQ and BEC models. In Section HI we explain the method of our 
analysis. Sections |IV| and [V| show the results of the analysis of mean multiplicities and particle 


spectra, correspondingly. Section VI concludes the paper. 


II. EXPLICIT TREATMENT OF HADRONIC GROUND STATES 


A. Ground state contribution 


Calculating particle multiplicities in a non-interacting quantum gas, one changes the sum¬ 
mation over discrete momentum levels by the integration over the continuous spectrum 
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Here gn is the degeneracy of the Pn-th momentum state, T is the system temperature, m is the 
mass of particles, and /i is the chemical potential. The integral over space coordinates gives 
the volume of the system V, and the factor (27r)3 appears since we use the natural units where 
h = 27Th and h = 1. For simplicity, we do not specify the physical character of the chemical 
potential p at the moment and consider only one type of particles. 


^ Note that we use the same acronym to denote both the Bose-Einstein condensation and our framework which 
explicitly includes the condensate formation. 
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If the chemical potential approaches the value of the particle mass, the integrand in the 
second line of Eq. Q stays constant for p —)■ 0, because the singularity of the denominator is 
cancelled by the integration measure d^p. However, the first term in the sum over quantum 
levels becomes inhnite in the limit /r —>■ m at p = 0 
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Therefore, at the onset of BEC, when /i —)■ m, one should keep the summation over the low 
momentum states and start the integration in Q at |p| > 0. One can show, however, that 
in the thermodynamic limit, H —)■ cxd, one may separate the p = 0 term only and start the 
integration from zero [l3], namely 
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The degeneracy of the zero level is the same as that of the continuous spectrum, go = g = 2s + l, 
where s is the particle spin. 


B. Chemical potentials in a hadron gas 

In the case of a hadron-resonance gas the above simple picture becomes more complicated — 
one should include: the sum over all existing hadron states, the integral over the mass spectrum 
of broad resonances, and the decays of resonances defined by their branching ratios. In addition, 
the chemical potential is different for each particle type i. In the case of the EQ model, the 
chemical potential becomes a linear combination of the particle’s electric, Qi, baryon, Hj, and 
strange. Si, charges 

= Qil^Q + BipB + Sips ■ (4) 

The chemical potentials pq, ps and ps are the same in the whole system. They control the 
conservation of electric charge, baryon number, and strangeness. In the midrapidity region at 
the LHC these quantities are negligible, so it is reasonable to assume that 

^ 0 . (5) 

The charge conservation is connected with the conservation of the difference of the quark and 
antiquark numbers. However, a fast expansion and cooling of the quark-gluon plasma may 
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lead to the effective conservation of the sum of qnarks and antiqnarks P]. It can happen if the 
hadronization is so fast, that the qnarks do not have enongh time to annihilate and lower their 
abnndances to the eqnilibrinm valnes at the chemical freeze-ont temperatnre. This leads to the 
appearance of additional non-eqnilibrinm chemical potential 

^NEQ ^ ^ 

where and are nnmbers of light {u, d) qnarks and antiqnarks, while iV* and Nl are the 
nnmbers of strange qnarks and antiqnarks in the ith hadron. For example, the non-eqnilibrinm 
chemical potentials for pions, protons, kaons, and A’s read: 


/^TT 5 3/iq , fJjK T ) /^A T ■ ('^) 


We neglect the contribntions from heavier qnarks, becanse they are snppressed by the Boltz¬ 
mann factor. 

The non-eqnilibrinm chemical potentials may be encoded in the fngacity factors of the form 


Ti = exp 


NEQ 

T 


/Q Is 


( 8 ) 


As a matter of fact, the non-eqnilibrinm hypothesis of hadron prodnction at freeze-ont was 
originally proposed in terms of Tj’s [1], therefore, we will follow the tradition and use jq and 
7 s instead of /ig and /i^. We have introduced the condensation term A^"cond given by Eq. (|^ in 
the latest version of SHARE [32] • The calculations done in this way are denoted below as the 
BEC case. 


III. FITTING STRATEGY 

A. Centrality selection 

There are several technical issues that we have to address now. Different particles are 
measured in different centrality bins. This leads to problems with the interpretation of the 
data. The authors of Ref. [13] solved this problem by redistributing heavy strange particles in 
smaller centrality bins, in which light particles such as pions are measured. However, we have 
found that this strategy introduces an uncertainty in the thermal analyses, which leads to large 
errors for the output parameters, especially for and q^. 

Therefore, herein we do the opposite and merge the particles measured with the hner cen¬ 
trality steps into larger centrality sets at which strange particles are measured. This gives us 
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six centrality intervals: 0-10%, 10-20%, 20-40%, 40-60%, 60-80%, and 80-90%. The corre¬ 
sponding errors are recalculated following the standard procedure of error propagation. 


B. Multiplicities vs. ratios 


The experiment provides both multiplicities (strictly speaking the rapidity densities) and 
particle ratios. Some authors prefer to £t particle ratios instead of multiplicities, in order to 
have fewer number of parameters and to get rid of the system’s volume. However, one should 
be careful with this procedure, especially for the non-equilibrium models, because with given 
n multiplicities one may construct n{n — l)/2 ratios. This gives quite substantial number of 
ratios for large n, which are correlated. 

In practice, one usually restricts oneself to a smaller and experimentally available set of 
ratios, and thus studies only a sub-domain of the whole system. This strategy affects the 
results, in particular for the NEQ models, where each particle has its own pre-factor Tj (|^. 
For example, the multiplicity of protons is approximately proportional to Tp = while that 
of pious to = 7 g. Therefore the ratio of protons to pions will give us the factor which is 
much less sensitive to the change of jg than the mean multiplicities of protons and pions. The 
ratio of lambda to kaon cancels the information on strangeness, Ta/T^ = (7g7s)/(7g7s) = 7g, 
and gives us as much information as proton to pion ratio, etc. Therefore, even if one takes all 
n(n — l)/2 ratios, one may get different temperatures and other parameters for the same set of 
n particles. 

Another important aspect of the analysis is that ALICE measures only the rapidity densities. 
Below, we argue that Eq. ([^ can be transformed in this case to the form 


AN 

Ay 


, ^norm T Acond; 
dy 


(9) 


where dV/dy is the volume of the produced particles per unit rapidity. Note that the volume 
does not cancel in the ratios of the rapidity densities in the case with the condensate included. 


C. Neutral pion effects on the fits 

The BEC model gives typically a very good £t for all observed particles. However, the 
calculated number of the not yet measured 7r° mesons is extremely large. This is so, because 
the £t suggests a high value of y^, at which the neutral pions, as the lightest particles, condense 
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first. This behaviour suggests that one should make a reasonable estimate of the vr^ multiplicity 
and use it as an input in the BEC £t. 

The isospin symmetry of strong interactions suggests that the number of vr'^’s must be the 
same as the number of positively or negatively charged pions. Although the isospin symmetry 
is in fact slightly broken (in particular, 7r° is lighter than vr^) we neglect this effect in order to 
get a lower bound for the amount of condensate, and use the constraint 


NAc) = - {NAc) + A-{c))- 


( 10 ) 


Equation (10) represents an estimate of a possible production rate. The error bars are recalcu¬ 
lated as follows 


ANAc) = ^ (AiV^+(c) + AiV^-(c)). (11) 

Since the experimental errors for the 7r~^ and 7t~ mean multiplicities are practically the same, 
we obtain almost the same error bars for 7r°. We note that it is much more difficult to measure 
7r°’s than charged pions. At present, the spectrum for > 700 MeV is known only. It gives 
just about 1/3 of all expected vr^’s. In the pt range where both A and A are measured, the 
latter has about twice bigger relative errors. The error bars for As also increase fast, while 
going from high to low px in the existing data m- Therefore, we think that ( [IT| ) is a reasonable 
estimate. 


IV. FIT OF HADRONIC MULTIPLICITIES 


We perform our analysis using the rapidity densities of A, p, p, K^, A, cj), and 
In order to include the effect of the neutral pions on the hts, we also include the estimate 
based on Eqs. (10) and ( pT| . This gives us 14 data points for the centralities 0-10%, 10-20%, 
20-40%, 40-60%, 60-80%, and 10 data points for the centrality bin 80-90%, where 5=*= and 
are not measured. We do not include the short-living A'* (892)° |39] — in this way we may 
check if it follows the pattern obtained for other long-living particles. By doing so we further 
examine the applicability of the single freeze-out concept, see our discussion in [32] . 

Our results are shown in Figs. Bi and[^ The volume determined from the hts is practically 
the same for the NEQ and BEC models, because the number of particles in the ground state 
appears to be relatively small. The EQ volume is larger, because T > 1 makes the system 
denser in NEQ and BEC, compared to EQ. Interestingly, the explicit treatment of the ground 
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FIG. 1: (Color online) The left figure shows the system volume, while the right figure shows the 
system temperature obtained as the fit to mean multiplicities of hadrons, which are produced in PbPb 
collisions at the beam energy ^snn =2.76 TeV. The BEC lines correspond to the non-equilibrium 
model with the Bose condensate in the ground state Q. The NEQ lines correspond to the non¬ 
equilibrium model without the Bose condensate, and the EQ lines are obtained for the equilibrium 
model. The NEQ lines are from the paper [32], while the EQ and BEC lines are from |35j . 

state lowers the gap in temperature between equilibrium and non-equilibrium models. These 
two versions of the model even coincide for the peripheral collisions, as it was shown before 
in [35]. 

The left panel of Fig. shows that the fif is much better for BEC than for EQ. Even 
the larger numbers of degrees of freedom in EQ (iVdof = 12 for c < 80% and iVdof = 8 for 
c = 80-90% in EQ, versus 10 and 6 in BEC, respectively) do not help to decrease the x^/^dof 
values. The agreement with the data for individual particles is also signihcantly better for BEC 
than for EQ. For example, the K*{892)^ resonance, which was not included in the £t, deviates 
from the data by more than 2cr in most central collisions in EQ, similarly as in Ref. 1271 , while 
it agrees with the data within la in BEC for each centrality. The successful fit of the short 
living iF*(892)° with the parameters obtained for long living particles is a strong argument in 
favor of the single-freeze-out approximation. 

The values of x^/^doi even drop below unity in BEC, which means that the experimental 
errors are probably too large. To check this behaviour, we calculated the y^/iVdof contours and 



















FIG. 2: (Color online) The non-equilibrium parameters jq on the left and 7 ^ on the right in BEG |35j 
and NEQ [35]. The corresponding EQ values are equal to unity by definition. 

found that the minimum of x^/^dof is very shallow in BEG. The right panel of Fig. shows 
the 10% deviation from the points with the best fit in the T — 7 ^ plane for BEG. The centrality 




FIG. 3; (Color online) The x^/^dof in BEG and EQ on the left, and the best fit values for 7 ^ and T 
in BEG on the right together with the error bars. The contours on the right figure indicate the area 
where x^/A^dof increases up to 10 %. 

increases from lower to higher temperatures, as in Fig. [T} One may see that the increase of 
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l^doi by 10% determines the area that covers a substantial part of the region marked by 
the error bars. This is so because there is a sudden rise in the values of x^/N^of for increasing 
7 g. This effect elongates the contours towards smaller values of One may also notice an 
anti-correlation between and T. 

The number of particles in the condensate in EQ, NEQ and BEG models can be calculated 
using Eq. (|^ and our results obtained with this formula are presented in Fig. The grey 
bands give the upper and lower bounds on the amount of condensate in BEG. The width of the 
bands is determined by the values of T and y^ taken along the contours shown in Fig. (other 
parameters are taken as the optimal). 




centrality [%] centrality [%] 


FIG. 4: (Color online) The ratio of neutral pions in the condensate on the ground state to the 

total number of neutral pions + .^norm is shown on the left figure. The same for charged 

pions is shown on the right figure. The BEG lines were obtained using Eq. while the EQ lines 


show the neglected contribution of the ground state in the EQ model, see Eq. (12). 


The condensate rate in the NEQ model is not shown, as it is very large. It means that 
using NEQ model, one should always take into account the ground state effect. Otherwise, the 
obtained y^ and y^ would be too high, while temperature too small. 

In BEG the condensate rate is 2% in the most central collisions and 7 — 8% in peripheral 
collisions (with a sudden increase to 19% in the most peripheral collisions). The plots for 7r° 


and TT^ are similar, because of the isospin symmetry that we demand with the use of Eqs. (10) 
and ( |TT| ). The upper bound is higher for in the most central collisions, because temperature 
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is lower there and the difference in masses of 7r° and becomes more important. 

The increase of the pion condensate contribution with centrality in the BEC model is un¬ 
expected, because jq and therefore decrease with centrality. However, it has been observed 
already before |13] that for smaller volumes the ratio of particles in the condensate increases if 
the system’s temperature is sufficiently high. Therefore, the decrease of volume may produce 
the increase in the condensate rate. 

We also note that the ratio of the 7r° and vr^ number in the zero momentum level (neglected 
in EQ calculations). 


N. 


EQ 


cond 


( 12 ) 


exp {m^/T^q) - 1 

to the total number of 7r° and tt^-s also increases for more than 10%. The number of particles 
in the ground state — 0.7 in (12) for all centralities, while decreases from about 

700 to 7 at the most peripheral collisions, see [26] for the vr^ data. Therefore the ratio rises 
from 0.1% to 10%. It means that the ground state should be included for the analysis of very 
peripheral heavy-ion or proton-proton collisions even in equilibrium models. 


V. SPECTRA 


A. Covariant Cooper-Prye formula with the condensate 


Equation ([^ can be formally obtained from ([^ by inserting the distribution function which 
contains a delta function component 

(2vr)3 


/(P) = 


exp 
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(13) 


in the integral over momentum ([^. 


In what follows we generalize Eq. (13) to the form suitable for hydrodynamic and kinetic- 


theory applications. At first we note that the number of particles N contained in the static 
box of volume V may be obtained as the integral over the phase-space distribution function 
f{x,p) = /(t,x,p), namely 

N = J (fxd^p f{x,p). (14) 

Hence, for the particles in a static box we may use the identification f{x,p) = {27i)~^f{p). 


If the particles are produced on the freeze-out hypersurface Eq. (14) is generalised to the 


11 












form 


N = / dS, 


f{x,p). 

-C/n 


(15) 


Then, the covariant generalisation of (13) for the systems which are boost-invariant and cylin- 
drically symmetric is 

1 9 


f{x,p) = 


(27r)3 exp (2^) - 1 


1 + ^^S{p-X)6{p-Y)S{p-Z) 


(16) 


exp rp 

Here V = / is the Lorentz invariant volume of the system, where the condensate is 

created. The quantity U is the four-velocity of the fluid element at freeze-out, while the four- 
vectors X, Y, and Z dehne the three spatial directions in the fluid local rest frame. In the 
center-of-mass system of the colliding nuclei, the form of the four-vectors X, Y , and Z is BSI: 


X = (sinh 6t cosh r/y , cosh 6t cos 0 , cosh 6t sin 0 , sinh 6t sinh r/y) , 
F = (0 , —sini/) ,cos(/) ,0) , 

Z = (sinhr^y ,0,0 ,coshr 7 y) . 


(17) 

(18) 
(19) 


Here Ot is the transverse rapidity of the fluid element, (j) is the azimuthal angle, and r/y is 
the space-time rapidity. In the local rest frame of the fluid element we have: X = (0,1,0,0), 
Y = (0,0,1, 0), Z = (0, 0, 0,1). The projections of the four-momentum along X, Y, and Z are 


px = p ■ X = rriT sinh 6t cosh(?/ — p\\) — Pt cosh 9t cos(0p — 0), 
Py = p-Y = Pt sin(0 - 0p), 

Pz = p ■ Z = m'rsinh(r 7 y — y). 


( 20 ) 

( 21 ) 

( 22 ) 


Here we used the parameterisation of the four-momentum: p^ = {p^, P^,p"^,P^) = 

{rriT cosh y, px cos 0p , px sin 0p , mx sinh y) where p^ = Ep = a/p^ -l- is the energy of a 
particle, y and (j)p are the particle rapidity and azimuthal angle of the momentum, px = 
is the transverse momentum, and mx = ^Jpx + ~ the transverse mass. 


The Lorentz covariant integral in (15) is 
1 


- - / 
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p (271)3 exp - 1 L 


1 + '^S(p-X)S(p-Y)S(p-Z) 


(23) 


Because of the Lorentz covariance and quadratic dependence of the Bose-Einstein distribution of 
three-momentum, one can conclude that N'^ = nU^ with the density n given by the expression 
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In the Cracow model, the hnal expression for AN (the number of particles in the rapidity 
interval Ay) is 

-^cond 


AN = / nU^ = Ay n„orm + 


— A f I ^cond 

= Ay— I nnorm “1“ 

dy \ V 
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Consequently, for the rapidity density one gets 


AN dV 

= ^ Ur, 


dV 


-N, 


cond‘ 
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(26) 


Ay dy Vdy 

In our case V = Tir'^^^TfAycond, where A|/cond is the rapidity range where the condensate may 
be formed. In this work we assume A?/cond = 1- Consequently we get 


AN _dV 

~~l ^norm “1“ tVcond- 

Ay dy 


(27) 


B. The particle spectra with the local condensate 


The momentum distribution is obtained from the formula 

„ dN dN 




^ d^p dy d(j)p pt dpx 
where the condensate part of the distribution function can be written as 

— d ( sinh 6t — 

m ■ 


SiPx) 5 {py) S{pz) = S (sinh^T - —) 5(0 - (j)p) 5(r?|| - y). 

Pt V m/ 


(28) 


(29) 


In our previous work we have found that the Cracow single freeze-out model can describe the 
LHC data for 2.76 TeV Pb-|-Pb collisions very well [3ll |32]. Therefore, at hrst we test effects 
of the Bose condensation using this model. In this case sinh(d'r) = ’"/'T/, where r and r/ are 
the freeze-out radius and proper time, respectively. Using this relation we may further write 

d{px) d{pY) S{pz) = - SlpT-r—] 6{p\\-y). (30) 

PTm \ Tf J 

Using this result we obtain the hnal expression for the momentum distribution 
dN - 
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where — r) is the Heaviside step function. 


The hrst term in the (31) is the usual one, while the second term with the delta function 
corresponds to the Bose condensate. The integrals over 0 and are cancelled by the delta 
functions and we obtain: 
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where A"cond is the number of particles in Bose condensate as in ([^. It follows from equation (32) 
that Bose condensate adds a constant number of particles to the usual momentum spectrum 
and these particles may have a momentum up to 


m 'x' max 

p™ = m - 


(33) 


We added the zero momentum level into THERMINATOR using the Dirac delta function 


(30). It turns out that in order to reproduce the spectra we may keep the same r-max/T/ ratio 
as in our previous paper |32] and just rescale r^ax and t/ in such a way that we obtain the 


volume determined from the studies of multiplicities described in Sec. |IV[ The results for pions 
are shown in Fig. The error bars for tt^ are not shown, because they are of the order of the 
symbol size. 

One may notice a very good overall agreement. In the BEG model the spectra of other 
particles including protons, kaons, K*{892)^ and 0(1020), also agree with data, and are almost 
identical to those that we obtained in j32]. The agreement with data for pions at high pt is 
even improved compared to [32j, because we use the same rf^a^/rf, but have bigger temperature 
in BEG than in NEQ. Therefore there are more high-p(r particles in BEG than in NEQ. 

On the other hand, the agreement of pions in the EQ model with the data is rather good, 
but other particles are described much worse [32j. For example, protons in EQ model deviate 
from data for more than the factor of two in the most central collisions at low pt [35]. 

The largest difference for pions is also seen at low pt- There is a step at pt = in 


the BEG model, see Eq. (33), because of the contribution of the particles from the ground 
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[GeV/c] [GeV/c] 

FIG. 5: (Color online) The pT spectra calculated in the BEC model (solid lines) and EQ model (dashed 
lines) at different centralities are compared with the existing data (symbols) for neutral pions |42j on 
the left and with the data for charged pions [2^ on the right. 


state. This is the artefact of the approximation that all particles from the low discrete levels 
are accumulated in the ground state ([^. The step is larger for neutral pions, because they 
are lighter and therefore condense earlier. Unfortunately, tt^’s have not been measured yet at 
low pt- There are only several points for charged pions, vr^, which are plotted in the linear 
scale without the 2” factors in Fig. One can see that the EQ lines do not describe the data 
in the low pt region for the very peripheral collisions. For the most central collisions the EQ 
lines go along lower edges of the error bars. This cannot be improved, because the EQ model 
overshoots protons at these centralities [32l [35] . Therefore it would be impossible to describe 
in EQ neither low px pions nor protons when the experiment provides the improved data with 
lower error bars. 

The BEC line with the step is within the error bars for the most central collisions, but the 
step is too big for more peripheral collisions. However, the smallest px points are htted at all 
centralities. One may also notice that the amount of extra particles to the left of the step is 
approximately equal to the number of missing particles to the right of the step. Therefore the 
inclusion of several more levels would lead to hner steps and smoother line, which goes within 
the error bars. 

The grey bands indicate the maximal and minimal amount of condensate at the 10% higher 
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FIG. 6: (Color online) The low momentum part of the spectra for charged pions in central and semi¬ 
central collisions on the left and in peripheral collisions on the right are compared to the calculations 
in the BEC (solid lines) and EQ (dashed lines) models. The grey bands correspond to minimal and 
maximal condensate rate shown in Eig. which is obtained assuming 10% increase of the y^/A^dof) 
see Eig. [fright. 


X^/iVdof from the best fit, see Fig. [fright and Fig. The grey band is within the error bars for 
most central collisions. One can even double the contribution from condensate there. However, 
in the most peripheral collisions the upper bound for the condensate is dehnitely too high and 
should be decreased by a factor of two^. Thus, the combined data on mean multiplicities and 
spectra are compatible even with the constant condensate rate at the level of 5%. 


VI. CONCLUSIONS 

In this work we have analysed the possibility of the pion condensation in heavy-ion collisions 
at the LHC energies. This phenomenon has attracted a lot of attention in the past (see [16] 
and references therein). 

Our study shows that the pion condensate rate in the present data is 2% in the most central 
collisions and 7 — 8% in peripheral collisions (with a sudden increase to 19% in the most 

^ One should remember that the discussed low-px region of tt^ spectra gives a small contribution to the mean 
multiplicities, therefore the spectra are much more sensitive to the condensate. 
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peripheral collisions). The obtained /^AoI valnes for the thermodynamic parameters have a 
shallow minimnm. It resnlts in a large nncertainty in the determination of the exact nnmber 
of pions in the condensate. The spectra show an opposite trend with centrality allowing larger 
condensate contributions in the most central collisions as compared to the peripheral ones. 

The data for charged pion spectrum at the 80-90% centrality would help to conhrm the 
increase of the condensate rate with centrality. However, the most important are the measure¬ 
ment of the 7r° multiplicities and the 7r° spectra in the low pt region. 
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